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Abstract 

The gravitational collapse of a star is an important issue both for general 
relativity and astrophysics, which is related to the well known "frozen star" 
paradox. This paradox has been discussed intensively and seems to have 
been solved in the comoving-like coordinates. However, to a real astrophys- 
ical observer within a finite time, this problem should be discussed in the 
point of view of the distant rest-observer, which is the main purpose of this 
paper. Following the seminal work of Oppenheimer and Snyder (1939), we 
present the exact solution for one or two dust shells collapsing towards a pre- 
existing black hole. We find that the metric of the inner region of the shell 
is time- dependent and the clock inside the shell becomes slower as the shell 
collapses towards the pre-existing black hole. This means the inner region of 
the shell is influenced by the property of the shell, which is contrary to the 
result in Newtonian theory. It does not contradict the Birkhoff's theorem, 
since in our case we cannot arbitrarily select the clock inside the shell in 
order to ensure the continuity of the metric. This result in principle may be 
tested experimentally if a beam of light travels across the shell, which will 
take a longer time than without the shell. It can be considered as the gener- 
alized Shapiro effect, because this effect is due to the mass outside, but not 
inside as the case of the standard Shapiro effect. We also found that in real 
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astrophysical settings matter can indeed cross a black hole's horizon accord- 
ing to the clock of an external observer and will not accumulate around the 
event horizon of a black hole, i.e., no "frozen star" is formed for an external 
observer as matter falls towards a black hole. Therefore, we predict that 
only gravitational wave radiation can be produced in the final stage of the 
merging process of two coalescing black holes. Our results also indicate that 
for the clock of an external observer, matter, after crossing the event hori- 
zon, will never arrive at the "singularity" (i.e. the exact center of the black 
hole), i.e., for all black holes with finite lifetimes their masses are distributed 
within their event horizons, rather than concentrated at their centers. We 
also present a worked-out example of the Hawking's area theorem. 

Key words: Black Holes, Classical Theories of Gravity, Spacetime 
Singularities 

PACS: 04.20.Jb, 04.70.Bw 



1. Introduction 

The "frozen star" is a well known novel phenomenon predicted by general 
relativity, i.e. a distant observer (O) sees a test particle falling towards a 
black hole moving slower and slower, becoming darker and darker, and is 
eventually frozen near the event horizon of the black hole. This process was 
vividly described and presented in many popular science writings [H, 0, 0] 
and textbooks jl, 0, @, 0, S B]- The time measured by an in-falling (with 
the test particle) observer (O') is finite when the test particle reaches the 
event horizon. However, for O it takes infinite time for the test particle 
to reach exactly the location of the event horizon. In other words, O will 
never see the test particle falling into the event horizon of the black hole. 
It is therefore legitimate to ask the question whether in real astrophysical 
settings black holes can ever been formed and grow with time, since all real 
astrophysical black holes can only be formed and grown by matter collapsing 
into a "singularity". 

Two possible answers have been proposed so far. The first one is that since 
O' indeed has observed the test particle falling through the event horizon, 
then in reality (for O') matter indeed has fallen into the black hole. It seems 
that this paradox could be solved in the point of view of O' or other comoving- 
like coordinates. However, since O has no way to communicate with O' once 
O' crosses the event horizon, O has no way to 'know' if the test particle 
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has fallen into the black hole. More importantly, in a real astrophysical 
observation within a finite time, it is more meaningful to discuss this problem 
in the coordinate of O, which is the main purpose of this paper. The second 
answer is to invoke quantum effects. It has been argued that quantum effects 



may eventually bring the matter into the black hole, as seen by O [10 



However, as pointed out by even in that case the black hole will still 
take an infinite time to form and the pre-Hawking radiation will be generated 
by the accumulated matter just outside the event horizon. (However, as we 
shall show in this work, the in-falling matter will not accumulate around the 
event horizon.) It has also been realized that, even if matter did accumulate 
just outside the event horizon, the time scale involved for the pre-Hawking 
radiation to consume all accumulated matter (and thus make the black hole 



"black" again) is far beyond the Hubble time 12] . Thus this does not answer 



the question in the real world. Apparently O cannot be satisfied with either 
answer. 

In desperation, O may take the attitude of "who cares?" When the test 
particle is sufficiently close to the event horizon, the redshift is so large that 
practically no signals from the test particle can be seen by O and apparently 
the test particle has no way of turning back, therefore the "frozen star" does 
appear "black" and is an infinitely deep "hole". For practical purposes O 
may still call it a "black hole", whose total mass is also increased by the 
in-falling matter. Apparently this is the view taken by most people in the 
astrophysical community and general public, as demonstrated in many well 
known textbooks jl, H, @, 0, S, @, E3] and popular science writings [H, 0, 0| . 
However when two such "frozen stars" merge together, strong electromag- 
netic radiations will be released, in sharp contrast to the merging of two 
genuine black holes (i.e., Schwarzschild singularities); the latter can only 
produce gravitational wave radiation Thus this also does not answer 

the question in the real world. As we shall show in this work this view should 
be changed, because in real astrophysical settings O (in the Schwarzschild 
coordinates) does observe matter falling into the black hole. 

This problem is closely linked to the problem of collapse. Great interests 
have been paid on the problem of gravitational collapse. Since the seminal 
work of Oppenheimer and Snyder (1939) [14|. numerous works have been 
done on this issue. For example, the result of [14] has been extended to more 
realistic and complicated cases, e.g. considering pressure, inhomogeneous, 



the effect of radiation and rotation, and so on (e.g. [15|, [l6|, UJ1, UM, Il9j ]) 



There are also a lot of works concerning collapsing shells. The paper of Isarel 



3 




Figure 1: The configuration of the dynamical case for one shell in the comoving coor- 
dinates, m and m s are the total gravitating masses of the black hole and the shell, 
respectively, a' and a are the radii of the inner and outer boundaries of the shell in the 
comoving coordinates, respectively. 



2o| found the equation of motion for one thin (without thickness) shell in the 
vacuum. Khakshournia and Mansouri [21] obtained the equation of motion 
for one shell in vacuum with thickness up to the first order. Gongalves 22 
discussed the solution for two thin shells in vacuum. However, in many real 
astrophysical settings, the matter shell collapsing towards a black hole should 
have a finite thickness and is much less massive than the black hole. In this 
paper we investigate such more realistic cases, and in particular to find the 
differences between these solutions and that in 14j, and the implication for 
the "frozen star" paradox. In Sec. [2J we present the dynamical solutions 
for one and two shells, in which the shell is spherically symmetrical and 
pressureless. We also show the evolution of event horizon (the boundary of 
the region from where the photon cannot escape to space-like infinity) and 
apparent horizon (the boundary of the trapped region, see the detail of the 



definition of event and apparent horizon in 13j and 23)) in the one shell 
case. In Sec. |3j we discuss the implications of these solutions and make our 
conclusions. We adopt G = c = 1 throughout this paper. 



2. The dynamical solution 

2.1. The solution for one shell 

Fig. [T] is the configuration of the one shell case. To investigate the 
dynamical solution, we follow the method in [lij ]. i.e. obtaining the solution 



4 



in the co moving coordinates first. The general form of the metric is 



ds 2 = dr 2 - e Q dR 2 - e^idQ 2 + sin 2 fl# 2 ), (1) 

where lj(R,t) and lj(R,t) are the unknown functions. In the comoving 
coordinates, the only non-zero component of the energy momentum tensor 
is T4 = p. As the results in we can obtain the following equations from 
the field equation 

e <" = (Ft + G) 4/3 , (2) 
e ° = e"w ,2 /4, (3) 

^P = \{t + G/F)-\t + G'/F>)-\ (4) 

where F and G are arbitrary functions of R. Dot and prime are the partial 
differentiation with respect to r and R, respectively. At r = 0, we obtain 

F = - v /67rp (G' 2 + C 1 ). (5) 



As in [14J, we choose G = i? 3//2 . 

In region I and III, po — 0, therefore, F = constant. We will determine 
this constant and C\ by the matching condition in the next paragraph. As 
obtained in pjj, in region III Fm = — §\/^o> where tq = 2M. According to 
the well-known Birkhoff's theorem for spherically symmetric mass distribu- 
tion, M here is the total gravitating mass of the system, i.e., including rest 
mass and gravitational energy. 

Throughout this paper, our matching conditions are chosen such that 
the metric is continuous across the boundaries between the different regions. 
In Sec. El we will verify that in our solutions found this way the extrinsic 
curvature is also continuous at the boundaries. Since Fm = — fy^To, to 
assure the metric is continuous at the boundary between region III and II, 
i.e. Fjn(R = a) = Fu(R = a), according to Eq. (jSJ), we obtain the value of 
C x as Ci = -^-M - a 3 . Therefore, to assure the metric is continuous at the 
boundary between region II and I, i.e. Fn(R = a') = Fi(R = a') 

of F in region I should be F\ = — |*/2[M — |7rp (a 3 - a' 3 )] = 



2 v L 3 

— fv^o ( we also denote m s = ^7rpo(a 3 — a' 3 )). In Sec. I2.1.2[ we will show 
that m and m s are the total gravitating masses of the central black hole and 
the shell, respectively. 

For the spherically symmetric ball case (just as the solution in [Hj]), C\ 
must be 0. However, in the case that a shell collapses towards a pre-existing 
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black hole, C\ can be any value. For example, G\ < means that the mass of 
the central black hole is smaller than the mass of the inner region interior to 
the shell in the spherically symmetric ball case. Thus, the gravitational force 
on the inner layers of the shell will become so weak that the inner layer will 
take more time than the outer layers (it can be seen from Eq. [7j) to arrive 
at r = [13]. Actually, it means the crossing of the layers will happen. We 
will not discuss the crossing case in this paper, since no comoving coordinates 
exist in that case. This means that we will only study the cases when G\ > 0. 
Therefore, our result cannot be compared directly with the result obtained 
for a collapsing shell (without a black hole in the center, i.e., C\ < 0) either 



without any thickness [20] or only with first-order thickness [2l[ . 

Up to now we have obtained the solution in the comoving coordinates. 
Then we try to transform the solution into the ordinary coordinates, in which 
the metric has the form as 

ds 2 = B(r, t)dt 2 - A(r, t)dr 2 - r 2 (d6 2 + sin 2 6d(f) 2 ). (6) 

The transformation of r is obvious. We must choose 

r = e" /2 = (Ft + G) 2/3 . (7) 

Using the contravariant form of the metric and requiring that the g tr term 
vanishes, we have 

t'/i = rr'. (8) 

Substituting the expressions of r in three regions into Eq. (jSJ), we can 
obtain three partial differential equations of t(R,r). 
In region III, the solution of Eq. (jSJ) is 

t = L(x), x = -^-{B?' 2 - r 3 / 2 ) - 2J¥F + r In ^ + ^ F \ ( 9 ) 



In region II, we note the solution of 

t'/i = -\^R [R 3/2 - T^poCR 3 + [v/6tt p {R 3 + d) - 67rp i? 3/2 r] 

(10) 

as t = M(y), where y is a function of R and r. 
The characteristic equation of Eq. (flOl) is 
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(11) 

In region I, the solution of Eq. (jSJ) is 

t = N(x'), x' = -^=(R 3/2 - r 3 / 2 ) - 2^7 Q + r' \n ^ + (12) 
3 V r o V^ _ V r o 

where L, M, and iV are arbitrary functions of their arguments. Since the 
metric in the ordinary coordinates in region III is the Schwarzschild metric, 
L(x) should be x. 

2.1.1. Analytic solution when C\ — 

Eq. (fTOj) can be solved analytically only when G\ = and the results are 

t = M(y),y = ^(i? 2 - a 2 ) + 1(1 - ^r)§, (13) 

where fc = 

Then if we require L = M at i? = a for any r, the form of M must be 
M(y) = ^=[a^- [kayf 2 ] - 2 + r In ^ + ^g , (14) 

where y = /(x') is determined by the following equation, 



x 



3-v/r, 



= [a'§ - (W) 3/2 (i/ - ^a' 2 + la 2 ) 3 / 2 ] - 2^a' [y - ±(a* - a 2 )] 

[y-±(a' 2 -a 2 )] + ^ 
Jka' [y - f(a' 2 - a 2 )] - ^ 



+r n In 



Using the relations 



1-r 2 ) 



and 



^ = t 2 (l-r^), 



.(15) 

(16) 
(17) 
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we can obtain the metric in the ordinary coordinates 
In region III, 

g a = l- r /r, (18) 
g rr = -(l-r /r)-\ (19) 



In region II, 



Since 



therefore, 



■(I-—)" 1 . (20) 



r (aky-r )(l - \Vkr)z 



(21) 



r (afcy-r ) 2 (l-f^r)3 kR 3 , 

9tt ~ ^F? ' (22) 

where k = §7rpo- 

In region I, in the ordinary coordinates, we can obtain 

fht = h(t)(l + C A /r), (23) 

g rr = -(l + C 4 /r)-\ (24) 

If we require g rr and g u are continuous at the inner boundary of the shell 
with Eq. (120!) and ff22|) . respectively, we have 

C 4 = -r' (25) 

and 

(1 - l^) 2 (UK 2 - * 2 ) + j(l - 1^)1] - a 2 } 2 
[(l-|v^r)l-r^/a'] 2 fc[±(a' 2 - a 2 ) + ±(1 - §V^r)i] 3 ' 

where r = r(t) is determined by the following equation 

2 [at - (fcay(r)) 3 / 2 ] - 2 y^FK + r In + v^ . ( 27 ) 



3 \/^ ' \fkay{r) - y/fZ 
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2.1.2. Numerical solution when C\ > 

In general cases, Eq. ([TO]) can be solved numerically only. From Eq. 
([9]), the value of t at the outer boundary of the shell can be determined. 
It can be easily shown that the value of t along every curve determined 
by Eq. flJTJ), i.e. r = t(R), is constant. Therefore, using the value of 
t from Eq. ([9]) as the boundary value, we integrate Eq. ( TTTj) along the 
characteristic lines and obtain how the shell evolves with coordinate time t. 
In Fig. [21 we show the evolution curves for the cases C% — and C% > 0. 
As it can be seen, the inner boundary of the shell crosses the Schwarzschild 

radius with finite time and approaches to a position r* = §rdr 3 — ~r' (This 
is obtained from the transformations Eq. [7] and Eq. |14j). whereas the 
outer boundary approaches to the Schwarzschild radius. (However, if there 
is another shell outside this shell, the outer boundary of this shell will also 
cross the Schwarzschild radius with finite time, as we will show in Sec. 12.2ft 
As shown in Fig. [2] the event horizon of the system increases just from 
the location r = 2m = r' Q . Therefore, r = r' Q corresponds to the event horizon 
of the central black hole and m is the total gravitating mass of the central 
black hole. As a result, m s can be recognized as the total gravitating mass 
of the shell. 

We find that the metric in region I is not the standard Schwarzschild 
metric; actually, it should be a function of time in the dynamical case. The 
evolution curve of h(t) is shown in Fig. [3j When the shell is far from the 
central black hole, the value of h(t) is almost equal to 1, i.e. the metric is ap- 
proximately the standard Schwarzschild solution. When the shell is near the 
horizon, the value of h(t) decreases rapidly and approaches to 0, i.e. the clock 
in region I becomes slower and slower. It does not contradict the Birkhoff's 
theorem, since in the proof of the Birkhoff's theorem the region outside the 
central mass is vacuum |4(. Therefore, in this case we could adjust the clock 
to ensure that metric asymptotically approaches to the Minkowski metric 
at infinity. However, in our problem there is a shell, but not the vacuum 
outside the region I. Thus we should ensure that the metric is continuous at 
the boundary between region I and the shell and therefore cannot arbitrarily 
select the clock inside the shell, i.e. region I. 

2.2. The solution for two shells 

In this section, we investigate the solution for two shells. The configura- 
tion of this problem is shown in Fig. [H The process to obtain the solution 
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Figure 2: The dynamical solutions for one shell case, (a) and (b) are evolution curves for 
a = 5ro, a' = 2.5ro, and Tq = l/8ro (i.e. C\ — 0) with comoving time and coordinate 
time, respectively. The evolution of the event and apparent horizons are also shown, (c) 
and (d) for a = 5r , a' = 2.5r , and r' — 0.5r (i.e. C\ > 0) with comoving time and 
coordinate time, respectively 
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r(R=a')/r 



Figure 3: The evolution of h(t) with the position of the outer boundary of the shell for 
the case a — 5ro, a' = 2.5ro, and r' Q = l/8ro. 




Figure 4: The configuration of the dynamical case for two shells in comoving coordinates. 
to, toi and TO2 are the total gravitating masses of the black hole, the shell 1 and 2, respec- 
tively. a[ and a\ are the radii of the inner and outer boundaries of shell 1, respectively. 
a' 2 and a2 are the radii of the inner and outer boundaries of shell 2, respectively. 
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is similar to the one shell case, however, there are five regions. The result of 
F is 



2 V^*0 5 ^ ^ ^1 



F = < 



-§V^(fl 3 -a?)+rg ;ai< J R<a 1 
"fv^o ; ai < < a' 2 



(28) 



where r = 2M, Tq = 2m, Tq = 2(m + mi). 

The form of the transformation in region V is the same as that in region 
III of the one shell case. The transformations for other regions are shown 
below. 

In region IV, 



t'/i = --VR [R 3 ' 2 - r v /6vrp 2 ( J R3 + C 2 )]-^ 6n p 2 (R 3 + C 2 ) - 67rp 2 R 3 / 2 r]. 

(29) 

The characteristic equation of Eq. (|29|) is 



% = \VR [R 3/2 - t^Qtt P2 (R3 + C 2 )]-i[ v / 67rp 2 ( J R3 + c 2 ) - Q-k P2 R 3 ' 2 
dri A 



where C 2 



4-rrp2 

In region III 



M-al. 



(30) 



t = 0(x"), x" 



3^o 



(R 3 / 2 - / 



3/2> 



2 v rr f Q + r In 



r + \/r r 



r — \/r P 



where 0(x") is an arbitrary function of x". 
In region II, 



(31) 



t'/i = -- v^R [i? 3/2 - r^6ir Pl (R 3 + d^—^^Gnp^R 3 + d) - 6tv Pi R 3/2 t}. 

(32) 

The characteristic equation of Eq. fl32|) is 

dr _ 2 
Iff ~ 3 

where C*i 



[R 3/2 - rv/evrp^/? 3 + W^ Pl {R 3 + d) - 67r PlJ R 3 / 2 r], 

(33) 

-^(m + m x ) - a 3 . 
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Figure 5: The dynamical solutions for two shells case, (a) and (b) are evolution curves 
for a 2 — 10r , a' 2 — 8r , a x = 5r , a[ — 2.5r , r' = l/3r , and j-q = 2/3r with comoving 
time and coordinate time, respectively. 



In region I, 

t = P(x'), x' = -^=(R 3 / 2 - r 3 / 2 ) - 2 y^7 + r' Q In ^ + , (34) 
3 V^o \A*-V r d 

where P(x') is an arbitrary function of a/. 

To obtain the evolution curves of the shells, we solve the above equations 
numerically with the same method in Sec. 12.1.21 Fig. [5] shows the evolution 
curves of two shells. The outer boundary of shell 1 crosses the Schwarzschild 
radius corresponding to the total gravitating mass within finite coordinate 
time. In Fig. [HI we show the evolution curves of the outer boundary of shell 
1 in the case with or without shell 2. As we have pointed out in Sec. 12.11 the 
clock inside the shell is slower, therefore the outer boundary of shell 1 will 
take more time to approach the asymptotic line in the case with shell 2. 



3. Discussions and Conclusions 

In this work, our matching conditions selection requires that the metric 
is continuous at the boundaries. Here we verify that in our solutions found 
this way the extrinsic curvature is also continuous at each boundary (see the 



detailed discussion about the extrinsic curvature in |23J]). Since the boundary 



is defined as R = constant, the normalized normal vector is 
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t/r„ 



Figure 6: The comparison of the evolution of the outer boundary of shell 1 between the 
case with (solid) or without (dashed) shell 2. The parameters of the shells are 02 = 10r , 
a' 2 = 6r , a[ = 2.5r , a\ = 5r , r' = l/5r , and j-q = 2/5r , where r is the Schwarzschild 
radius corresponding to the total gravitating mass of the system, i.e., the sum of the 
gravitating masses of the pre-existing black hole, shell 1 and shell 2. The inset is the 
ratio of the solid line to the dashed line. It is seen clearly that outside matter, though 
spherically symmetric, does influence the motion of matter inside. 



n « = (°>^17^ > )- ( 35 ) 

The definition of extrinsic curvature is K a b = n Q;/ 3e"ef. We choose y a = 
(t, 9, ip) as the coordinates on the boundary. Therefore, the non-vanishing 
components of K a b are 

K TT = n T]T = -r*n R = ^n R g RR g TT , R = 0, (36) 
K ee = n . e = -T% e n R = ^g RR gee, R n R = (Ft + G) 2/3 , (37) 



K w = = -T R v n R = \g RR g w ^ R n R = (Ft + G) 2/3 sin 2 6. (38) 

Obviously, they are all continuous at all boundaries. 

With the dynamical solutions presented in this work, we conclude that 
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(a) In general relativity, the outer mass influences the inner space-time 
even if the mass distribution is spherically symmetric; 

(b) The metric inside the shell can not be the standard Schwarzschild 
form and is time dependent in the dynamical case; 

(c) A dust shell of finite thickness can indeed cross the black hole's event 
horizon, as observed by an external observer, except for its outer surface. 
This means that asymptotically only an infinitesimal amount of matter re- 
mains just outside the event horizon. Therefore such a "star" is qualitatively 
different from the so-called the "frozen star", which supposedly holds all 
in-falling matter outside its event horizon; 

(d) In the two-shell case, even the outer boundary of the inner shell can be 
observed by an external observer to collapse into the Schwarzschild radius, 
thus no matter can be accumulated outside a black hole's event horizon. 
Since we can mimic the outer shell as all matter between the observer and 
the in-falling shell being observed, we can conclude that the outside observer 
indeed "sees" the inner shell falls into the black hole completely. We note 
that the outer layer of matter does not necessarily block the view of the 
outside observer to the inner layer, since the outer layer can be optically thin 
to at least a certain wavelength of electromagnetic radiation, or is made of 
dark matter and is thus completely transparent to electromagnetic radiation. 
Even if the outer layer is optically thick completely, such as the case when 
the fall-back matter approaches to the black hole after the core collapse 
of a massive star, neutrinos can still find their way out and photons may 
eventually escape after a lengthy radiative transfer process. 

In our solution, the metric between the shell and the black hole should 
not be the standard Schwarzschild form, and there should be an additional 
time-dependent factor h(t) in g tt (please refer to equations Eq. (1251) and 
( I26p . Fig. [31 This factor h(t) decreases from 1 to when the shell collapses 
towards the Schwarzschild radius. This means that the clock in the inner 
region of the shell is slower compared with the case without the shell, and 
becoming even slower when the shell is falling in. This effect may be testable 
experimentally in principle, e.g., if a beam of light travels across the shell and 
a parallel beam of light passes outside the shell (far enough from the shell), 
then the two beams of light will travel at different velocities with respect to 
the observers outside the shell. This previously unrecognized effect is caused 
by the outer mass's influence to the inner space-time; such effect does not 
exist in the Newtonian Gravity. We could recognize it as the generalized 
Shapiro effect, because this effect is due to the mass outside, but not inside 
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as the case of the standard Shapiro effect. We comment in passing that 



in 



22], the author directly identified the metric between two thin shells as 
Schwarzschild metric and therefore his result can recover to our one shell 
solution only when the mass of the outer shell vanishes. 

In the one shell case, the shell could almost cross the Schwarzschild radius 
except the outer boundary; such a star is already qualitatively different from 
a "frozen star". Nevertheless, a shell can indeed cross the Schwarzschild 
radius in the two shell case, and the external observer can indeed observe 
a complete shell falling into a black hole, since the outer shell in general 
cannot block the view to the inner shell completely, as discussed above. In 
this sense we could observe the matter falls into a black hole and the "frozen 
star" paradox discussed in 10] is naturally solved. As pointed out previously 
by us, the origin of the "frozen star" is the "test particle" assumption for 
the in-falling matter, which neglects the influence to the metric by the in- 
falling matter itself 24|. In fact, for some practical astrophysical settings, 
the time taken for matter falling into a black hole is quite short, even for 



the external observer 24]. Therefore with full and self-consistent general 
relativity calculation, we find that the in-falling matter will not accumulate 
outside the event horizon, and thus the quantum radiation and Gamma Ray 
bursts predicted in and [12j are not likely to be generated. We predict 
that only gravitational wave radiation can be produced in the final stage 
of the merging process of two coalescing black holes. Future simultaneous 
observations by X-ray telescopes and gravitational wave telescopes shall be 
able to verify our prediction. 

It is also interesting to note that, as can be seen from Fig. [2b, [2H and 
[5b, in ordinary coordinates, the matter will not collapse to the singularity 
(r = 0) even with infinite coordinate time (if r' = 0, the inner boundary 
will take infinite time to arrive at r = 0). It means that in real astrophysics 
sense, matter can never arrive at the singularity (i.e., the exact center of the 
black hole) with respect to the clock of the external observer. Therefore, no 
gravitational singularity exists physically, even within the framework of the 
classical general relativity. 

As predicted by Hawking's area theorem, the event horizon grows monotonously 
and continuously from r' to ro when the shell collapses towards the black hole, 
whereas the apparent horizon jumps to the Schwarzschild radius when all the 
mass of the shell crosses into the Schwarzschild radius (see Fig. [2^, and 4b). 
Therefore, we have presented a worked-out example of the area theorem. 
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